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We present the results of a deep spectral analysis of all Swift observations of Mrk 421 between 
April 2006 and July 2006, when it reached its highest X-ray flux recorded until the end of 2006. 
We completed this data set with other historical X-ray observations. We used the full data set to 
investigate the correlation between the spectral parameters. 

We found a signature of stochastic acceleration in the anticorrelation between the peak energy 
(E p ) of the spectral energy distribution (SED) and the spectral curvature parameter (b). We 
found signature of energetic budget of the jet in the correlation between the peak flux of the SED 
(S „) and E p . Moreover, using simultaneous Swift UVOT/XRT/BAT data, we demonstrated, that 
during the strongest flares, the UV-to-X-ray emission from Mrk 421 requires that the curved 
electron distribution develops a low energy power-law tail. 

The observed spectral curvature and its anticorrelation with E p is consistent with both stochastic 
acceleration or energy-dependent acceleration probability mechanisms, whereas the power-law 
slope of XRT-UVOT data is close to that inferred from the GRBs X-ray afterglow and in 
agreement with the universal first-order relativistic shock acceleration models. This scenario 
implies that magnetic turbulence may play a twofold role: spatial diffusion relevant to the first 
order process and momentum diffusion relevant to the second order process. 



X 



The Extreme sky: Sampling the Universe above 10 keV - extreme sky 2009, 
October 13-17, 2009 
Otranto (Lecce) Italy 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http: //pos. sissa. it/ 



The Large Flares ofMrk 421 in 2006 



Andrea Tramacere 




Ep (keV) 

Figure 1: Scatter plot of the curvature (b) vs. E p . Solid red circles represent data from Tramacere et al. 
(2007). Black boxes represent Swift data from Tramacere et al. (2009), without the cooling-dominated 
events. Empty circles represent the whole XRT data set presented in Tramacere et al. (2009). 

1. Introduction 

BL Lac objects are Active Galactic Nuclei (AGNs) characterized by a polarised and highly 
variable nonthermal continuum emission extending from radio to y-rays. In the most accepted 
scenario, this radiation is produced within a relativistic jet that originates in the central engine 
and points close to our line of sight. The relativistic outflow has a typical bulk Lorentz factor 
of T « 10, hence the emitted fluxes, observed at an angle 9, are affected by a beaming factor 
6= 1 /(IX 1-/3 cos (9)). 

The Spectral Energy Distribution (SED) of these objects has a typical two-bump shape. According 
to current models, the lower-frequency bump is interpreted as synchrotron (S) emission from highly 
relativistic electrons with Lorentz factors y in excess of 10 . This component peaks at frequencies 
ranging from the IR to the X-ray band. In the framework of the Synchrotron Self Compton (SSC) 
emission mechanism, the higher-frequency bump can be attributed to inverse Compton scattering 
of synchrotron photons by the same population of relativistic electrons that produce the synchrotron 
emission (Jones et al., 1974). 

Mrk 421 is classified as a High energy peaked BL Lac (HBL) (Padovani & Giommi, 1995) be- 
cause its synchrotron emission peak ranges from a fraction of a keV to several keV. With its red- 
shift z - 0.031, it is among the closest and most well studied HBLs. In spring/summer 2006, 
Mrk 421 reached its highest X-ray flux recorded until that time. The peak flux was about 85 
milli-Crab in the 2.0-10.0 keV band, and corresponded to a peak energy of the spectral energy 
distribution (SED) that was often at energies higher than 10 keV. 

In this paper we use Swift UV and X-ray data (Tramacere et al., 2009) of the 2006 flaring activity, 
completed with historical X-ray observations, to interpret the correlation between the spectral pa- 
rameters in terms of acceleration processes and energetic budget of the jet. We remand the reader 
to Tramacere et al. (2009) that paper for a more complete picture. 
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Swift Xrt Apr-Jill 2006 

Tramacere et al. 2007b (data from 1997 to 2005) 





Figure 2: Scatter plot of L p vs E p . Solid red circles represent data from Tramacere et al. (2007) spanning 
from 1997 to 2005. Solid blue boxes represent Swift data from Tramacere et al. (2009). The green dashed 
dotted line represents the best fit using a power-law and the black dashed line represents the best fit using a 
broken power law. 

2. Spectral parameter trends 

Mrk 421 showed in the past several major flaring episodes marked by flux variations that went 
along with significant spectral variations (Massaro et al., 2004). The X-ray spectral shape in gen- 
eral exhibited a marked curvature that is described well by a log-parabolic model (Massaro et al., 
2004; Tramacere et al., 2007), that can be expressed in terms of the SED peak energy (E p ), of the 
SED peak flux (S p ), and of the peak spectral curvature (b) as: S(E) = S P I0~ b ( 1 °s(e/e„)) 2 _ 
In the following we investigate the signature of stochastic acceleration in the anticorrelation be- 
tween E p and b. Moreover, we look for signatures of energetic budget of the jet in the correlation 
between S p and E p . 

2.1 E p - b trend 

In Fig. 1 we show the E p - b trend for both XRT 2006, and historical data. Thanks to the 
very high statistics exhibited during the 2006 flaring activity, we were able to perform an highly 
temporally-resolved spectral analysis and to identify cooling-dominated observations (empty cir- 
cles in Fig. 1), characterized by a strong flux decrease and strong spectral softening. These cooling- 
dominated states are disentangled from the full dataset because they can bias the acceleration- 
driven trend. 

The E p -b scatter plot shows a clear anticorrelation, with the peak energy of the S component 
increasing as its spectral curvature is decreasing. This trend hints for an acceleration processes 
producing curved electron distributions, where the curvature decreases as the acceleration becomes 
more efficient. A first possible scenario is that in which the acceleration probability of the particle 
is a decreasing function of its energy (Massaro et al., 2004). An alternative explanation is provided 
by the stochastic acceleration (SA) framework (Kardashev, 1962; Tramacere et al., 2007). In the 
SA mechanism the momentum-diffusion term (D) is responsible for the broadening of the electron 
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energy distribution (n(y)), hence it is inversely proportional to the curvature in the observed X-ray 
photons (b) and to the curvature in the distribution of the emitting particles (r) . According to the 
relation between r and b (Massaro et al., 2004, r ~ 5b) we derive the following relation among the 
peak energy of the electron distribution (y p ), b, and E p : 

log(£ p ) = 2 1og( rp ) + 3/(5£). (2.1) 

Clearly, this trend predicts the anticorrelation between E p and b, in agreement with the observa- 
tional data reported in Fig. 1. 

2.2 S p -E p trend 

The Sp-Ep trend (Fig. 2 ) demonstrates the connection between the average energy of the 
particle distribution and the power output of the source. To obtain a deeper understanding of the jet 
energetics, we plot on the y axis L p = S p 4nD 2 L , where Dl - 134 Mpc is the luminosity distance 1 . 
In the case of synchrotron emission, we expect: S p « n(y^ p )B 2 8 4 and E p <x y\ p B6, where y^ p is the 
peak of n{y)y i , B is the magnetic field, and 5 is the beaming factor. 

It follows that the dependence of S p on E p can be expressed in the form of a power-law: S p <xEp a ' . 
The simple power-law fit gives a value of a - 0.42 ± 0.06, this value clearly rules out as main drivers 
both B(a = 2) and 8{a - 4), indicating y p as the main driver. 

A more detailed analysis of the scatter plot reported in Fig.2 shows that the trend has a break at 
about 1 keV. The broken power law fit gives two slopes of a\ — 1.1 ±0.2 and = 0.27 + 0.07, 
respectively. This break in the trend implies that for E p < 1 keV and L p < 10 45 erg/s, the driver 
follows the relation with a =* 1.0 (we define this state the quiescent sate), whilst for E p > 1 keV and 
L p > 10 45 erg/s, the driver relates to a - 0.2 (we define this to be the high state). 
A possible interpretation is that the break of S p - E p depends on the modulation of the number 
of emitted particles (N = Jn(y)dy). Since B and 8 have been ruled out as main drivers of the 
Sp-Ep trend, we assume that they have a small variance. In this scenario, as yj p (namely E p ) is 
increasing N is constant up to a maximum value of the electron energy density u e = Jym e c 2 n(y)dy 
corresponding to a L p of - 10 45 erg/s, that could represent the maximum energetic budget of the 
jet. Above the value of E p - 1 keV, N decreases to don't exceed the maximum energy content. 

3. UV-to-X-ray spectrum and SED modelling 

Combined simultaneous UVOT and XRT observations (Tramacere et al., 2009) shows SEDs 
(see left panel of Fig. 3 ) that can be classified in three categories: a) described by a log-parabola, 
b) described by a power law, c) described by a spectral law that is a power law at its low energy 
tail, becoming a log-parabola function at its high energy one (LPPL) 

The power-law spectral index, in the vF(v) representation, for model b and c were found in the 
range a v - [0.25 -0.4]. According to the standard synchrotron theory, these values are not com- 
patible with a high value of the minimum energy of the electron distribution (y m! „). Indeed, in this 
case we would expect a v - 4/3, a value much harder than the observed one. This implies a low 
energy power-law tail in n{y). Using the well known relation (Rybicki & Lightman, 1979) between 

1 We used a flat cosmology model with Hq = 0.71 km/(s Mpc) CIm = 0.27 and fl^ = 0.73. 
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Figure 3: Left panel: Three different spectral shapes of the UV-to-X-ray data from Tramacere et al. (2009). 
Red boxes represent a PL spectrum observed on 2006 June 16. Blue diamonds represent a LPPL spectrum 
on 2006 April 25. Green circles represent a LP spectrum on 2006 June 14. Right Panel: SSC model. Solid 
red circles represent data from, 06/23/2006 UVOT-XRT-BAT observations. Green up triangles shows UVOT- 
XRT data on 03/31/2005 from Tramacere et al. (2007). Empty grey polygons represent non-simultaneous 
EBL corrected TeV data (Albert et al., 2007; Yadav et al., 2007). The solid red down triangle represent a 
Whipple observation on 18,19 and 21 June 2006 (Lichti et al., 2008). Best fit parameters for the 06/23/2006 
state: source size R = 2.1 x 10 15 cm, 5=0.1 G,6 = 25.0, N=15.0 cm" 3 , s=23, r=0.65, y c .=2.85xl0 5 . 



the electron energy distribution spectral index (s ) and that in the S photons (s = 3 - 2a v ) we obtain 
s - [2.2 - 2.5]. A phenomenological description of n(y), consistent with the observed UV-to-X-ray 
spectrum, is given by the following spectral law: 

n(y) = K{yly c y s , 7<7c 

n(y) = K (y/y c T (s+r lo ^»,y > y c , (3. 1) 

where r is the curvature and y c the turn-over energy. The presence of a power-law feature and the 
range of observed spectral indices are relevant both in the context of Fermi first-order acceleration 
models and from an observational point of view. Indeed, our range of spectral index values (s - 
[2.2 - 2.4] ) is consistent with the predictions of the relativistic-shock acceleration (s - 2.3), both 
in the case of analytical or numerical approaches (Achterberg et al., 2001; Blasi & Vietri, 2005). 
Moreover, we note that Waxman (1997), studying the the afterglow X-ray emission of y-ray bursts 
(GRB), inferred an electron distribution index of s - 2.3 ±0.1. Although the emission scenarios 
assumed in the case of GRBs are more complex w.r.t. the one-zone homogeneous SSC model, the 
similar values of the electron energy spectral index hints for the relativistic-shock acceleration as an 
underlying acceleration mechanism common both to GRBs and HBLs. We remind that the power- 
law feature can also be consistent with a purely stochastic scenario, but in this case, it requires a 
fine tuning of the ratio of the acceleration to the loss timescale (s - 1 + t acc /t esc ) to get a universal 
index. 

An example of SSC modelling obtained using a LPPL electron distribution is reported in the right 
panel of Fig. 3 (for a more detailed analysis of the SED modelling we remand to Tramacere et al. 
(2009)). The red points represent the simultaneous Swift data on 06/23/2006, and the green points 
data for the 31/05/2005 state. Here we just stress the huge difference in the spectral shape and flux 
levels between the two states. Indeed, for the 2005 data the S component is well described both in 
the UV and X-ray window, by using as n(y) a pure log-parabola. On the contrary, the 2006 state 
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requires a LPPL electron distribution to fit the UV-to-X-ray data. Because of the dominant galaxy 
contribution below UV energies, 2005 Swift data don't allow to understand whether or not a low 
energy PL feature develops. We remark that such a feature could be tested comparing the prediction 
from a one zone homogeneous SSC model with simultaneous Swift and Fermi-LAJ MeV/GeV 
data. 

4. Conclusion and Discussion 

We have shown that the X-ray flaring activity of Mrk 42 1 results in a complex spectral evo- 
lutions due to drastic changes in the electron energy distribution probably related to a complex 
acceleration scenario. Indeed, in our analysis we found both signatures of first and second order 
acceleration processes acting at the same time. The E p -b trend is consistent with a SA scenario 
with the X-ray spectral curvature related to the acceleration rather than to the cooling process. 
The presence of a power-law low-energy tail, found during the brightest X-ray flares in 2006, and 
the corresponding values of the electron distribution index (s =* [2.2 - 2.4]) are consistent with the 
predictions of relativistic Fermi first-order acceleration models (s 2.3). Our findings hint for a 
simultaneous role of the first and second order processes both related to the magnetic field tur- 
bulence. The stochastic acceleration hence is related to a momentum diffusion coefficient which 
drives the curvature consistently with the E p -b observed trend. The first order acceleration, ob- 
servationally signed by the low-energy power-law spectral index, is linked to the spatial diffusion 
coefficient. Interestingly, recent particle-in-cell simulations by Spitkovsky (2008) obtain electron 
distributions that are compatible with this scenario. 

The S p - E p trend demonstrates the connection between the average energy of the particle dis- 
tribution and the power output of the source. The observed values of the expected power-law 
dependence (S p « E a p ) exclude B, and 8, indicating y p as the main driver of the S p -E p trend. The 
break in the S p -E p scatter plot (see Fig. 2) at about 1 keV, where the typical source luminosity is 
about hp - 10 45 erg/s, can be interpreted as an indicator of the energetic content of the jet. 
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